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ABSTRACT 

With orbital periods of the order of tens of minutes or less, the AM Canum Venatico- 
rum stars are ultracompact, hydrogen deficient binaries with the shortest periods of 
any binary subclass, and are expected to be among the strongest gravitational wave 
sources in the sky. To date, the only known eclipsing source of this type is the P = 28 
min binary SDSS J0926+3624. We present multiband, high time resolution light curves 
of this system, collected with WHT /ULTRAC AM in 2006 and 2009. We supplement 
these data with additional observations made with LT/RISE, XMM-Newton and the 
Catalina Real-Time Transient Survey. From light curve models we determine the mass 
ratio to be q = M 2 /M 1 = 0.041 ± 0.002 and the inclination to be 82.6 ± 0.3 deg. We 
calculate the mass of the primary white dwarf to be 0.85 ± O.O4M and the donor 
to be 0.035 ± 0.003M Q , implying a partially degenerate state for this component. We 
observe superhump variations that are characteristic of an elliptical, precessing ac- 
cretion disc. Our determination of the superhump period excess is in agreement with 
the established relationship between this parameter and the mass ratio, and is the 
most precise calibration of this relationship at low q. We also observe a quasi-periodic 
oscillation in the 2006 data, and we examine the outbursting behaviour of the system 
over a 4.5 year period. 

Key words: stars: individual: SDSS J0926+3624 — stars: binaries : close — stars: 
white dwarfs — stars: cataclysmic variables 



1 INTRODUCTION 

The AM Canum Venaticorum (AM C Vn) stars are ultra - 
compact binaries with periods from 5.4 (|Roelofs et alj|20ich 
to 65 m inutes and opti c al spectra dominat ed by helium 
(see, e.g. iNelemansl 120051 : iRamsav et al .1 [20071 for recent re- 
views). These systems consist of a white dwarf accreting 
matter via a helium accretion disc from a significantly less 
massive and hydrogen deficient donor star. In order to fit 
within the Roche lobe it is necessary for this donor to 
also be at least partially degenerate. AM CVn stars of- 



fer new insights into the formation and evolution of bi- 
nary star systems, with the short periods implying at 
least one common envelope phase in the history of the bi- 
nary, and the chemical composition suggesting helium white 
dwarfs o r CVs with evolved second aries as possible pro- 
genitors jN clcmans et al. 2001, 20icl. see also iMarsh et all 
2010; iKulkarni fc van Kerkwiikl 120101 for the recent identi- 
fication of a possible AM CVn progenitor). Close double- 
degenerate binaries are also one of t he proposed progenitor 
popu l ations of Type la supernovae l|Tutukov fe Yungelsonl 
Il98ll : IWebbink|[l984l : llben fc Tutukovlll984l ) and sublumi- 
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nous events (|Perets et al.ll2ofot ) . Finally, the mass transfer in 
these systems is thought to be driven by angular momentum 
loss as a result of gravitational radiation. Due to their very 
short periods they are predicted to be among the stronges t 
gravitational wave sources in the sky (|Nelemans et al.f 2004). 
and are the only class of binary with examples already 
known which will b e detectable by the grav i tational wave 
obser vatory LISA (|Stroeer fc Vecchiol 120061 ; iRoelofs et ail 
l2007h . 

Gravitational radiation has a huge influence on AM 
CVn systems, driving their evolution and determining their 
orbital period distribution, luminosities and numbers. De- 
generate stars expand upon mass loss and so stable mass 
transfer via Roche lobe overflow causes an evolution towards 
longer periods. The combination of decreasing donor mass 
and lengthening orbital period leads to a rapid decrease in 
the magnitude of the gravitational wave losses over time. 
There is therefore a significant drop in the mass transfer 
rate over the observed perio d range of the AM CVn popula- 
tion l|Nelemans et al.|[200lh . If the donor stars in AM CVn 
stars were completely degenerate then their masses would 
be a unique function of orbital period, and their mass trans- 
fer rates a function of the accretor mass and orbital period. 
However, the three current paradig ms for the binary for- 
mation path (w hite dwarf mergers, Nelemans et al.l l200ll ; 



ex-heli um stars, Iben & Tutukov 1991; CVs with evolved 
donors, Podsiadlowski et al. 2003) all imply partial degener- 
acy, to different degrees. A partially degenerate star must be 
more massive than a degenerate star to fit within a Roche 
lobe at a given orbital period. A less degenerate donor there- 
fore implies higher gravitational wave losses and a higher 
mass transfer rate. A test of the degeneracy of the donor star 
requires accurate mass determinations which have proved 
elusive, although some constraints were obtained for five 
systems using para llax measurements obtained with HST 
|Roelofs et al.ll2007l ). 

The p rototype AM CVn system was discovered 40 
years ago ( Sm all Tl967l : iPaczvhskl 1 19671 ) . but to date only 
~25 furt her objects of this class have been discovered 
(see, e.g., IRoelofs et all 120051; lAnderson et all 1 20051 . 120081; 



IRoelofs et al.H2009l ; iRau et all 120101 ) . One of these was the 
eclipsing system SDSS J0926+3624 l|Anderson et al.ll2005l . 
SDSS 0926 hereafter). SDSS 0926 is currently the only 
eclipsing AM CVn known, and has a period of 28 min, with 
eclipses lasting ~1 min. The mean fl-ba nd magnitude of this 
system is ~19.3 (I Anderson et al.ll20 05). but there is consid- 
erable out-of-eclipse variation, characteristic of the super- 
humping behaviour seen in many AM CVns and CVs which 
is attributed to the precession of an elliptical accretion disc 
|Whitehurstlll988! ; lLubowlll99ll ; [Smipson fc Woodlll998l ). 

In 2006 and 2009 we took high time resolution obser- 
vations of SDSS 0926 with the fast CCD camera ULTRA- 
CAM. The aim of these observations was to determine pre- 
cise system parameters for this system, using techniques we 
have in the past s uccessfully applied t o normal CVs (e.g 



nave m tne past successfully appl ied t o normal UVs (e.g. . 
Feline et al.l 120041; iLittlefair et al.ll2008l; iPvrzas et al.ll2009l : 



Southworth et all 120091 ; Copperwheat et alT 20ld ). Precise 
masses enable us to determine the degree of degeneracy of 
the donor star, and eclipse timings can be used to determine 
the angular momentum losses. We present these photomet- 
ric observations in this paper, as well as additional data 



collected with the Liverpool Telescope, XMM-Newton and 
the Catalina Sky Survey. 



2 OBSERVATIONS 
2.1 WHT/ULTRACAM 

In 2006 and 2009 observations of SDSS 0926 w ere made with 
the h igh speed CCD camera ULTRACAM |Dhillon et all 
2007) mounted on the 4.2m William Herschel Telescope 
(WHT). The 2006 observations were mainly taken over a 
three day period in the beginning of March. Weather con- 
ditions were reasonable, with seeing ~1" and good trans- 
parency. The 2009 observations were taken over three nights 
in January, and conditions for these winter observations were 
on the whole poorer, with variable seeing and transparency. 
Due to conditions, only a small number of orbital cycles 
were observed on two of the three nights. ULTRACAM is 
a triple beam camera and all observations were made using 
the SDSS u' , g' and r' filters. Average exposure times were 
~3s in 2006 and 1.8s in 2009. The longer exposure time was 
necessary for the 2006 observations due to the relatively low 
S/N of the tt'-band data. Our 2009 observations took advan- 
tage of a new feature in the ULTRACAM software, in which 
multiple u'-band exposures can be coadded on the CCD be- 
fore readout. Two coadds were used for the majority of the 
data, giving a w'-band exposure time of 3.6s, although this 
was increased during poor conditions. The dead time be- 
tween exposures for ULTRACAM is ~25ms. The CCD was 
windowed in order to achieve this exposure time. A 2 x 2 
binning was used in most of the 2006 observations to com- 
pensate for conditions. A complete log of the observations 
is given in Table [1] 

All of these data were reduced with aperture photom- 
etry using the ULTRACAM pipeline software, with de- 
biassing, flatfielding and sky background subtraction per- 
formed in the standard way. The source flux was deter- 
mined using a variable aperture (whereby the radius of the 
aperture is scaled according to the FWHM). Variations in 
transparency were accounted for by dividing the source light 
curve by the light curve of a nearby comparision star. The 
stability of this comparison star was checked against other 
stars in the field, and no variations were seen. We determined 
atmospheric absorption coefficients in the u' , g' and r' bands 
and subsequently determined the absolute flu x of our tar- 
gets using observations of standard stars (from lSmith et all 
2002) taken in evening twilight. We used this calibration for 
our determination of the apparent magnitudes of the source, 
although we present all light c urves in flux units determined 
using the conversion given in lSmith et al.l (|2002l ). Using our 
absorption coefficients, we extrapolate all fluxes to an air- 
mass of 0. For all data we convert the MJD (UTC) times 
to the barycentric dynamical timescale, correcting for light 
travel times. 



2.2 LT/RISE 

We supplemented our WHT/ULTRACAM data with addi- 
tional observations taken through the first half of 2009 with 



the high speed RISE camera (jSteele et al.l 2008) mounted 
on the Liverpool Telescope. These observations were taken 
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Table 1. Log of the observations. 
Night UT Exposure Number of 





start 


end 




time (s) 




Binning 


orbits 


Comments 
















WHT/ULTRACAM 




1 c:+ Mar 90Dfi 


22 


28 


04 : 


48 


3-4 




2x2 


13 


Qt_-Llllli _1_ , SUIIIL. jJajLUIlL-o Ul LIU UU 


2nd Mar 


20 


04 


04 : 


49 


2-3 




1x1/2x2 


18 


Seeing ~ 0.8", clear. 


3rd Mar 


19 


56 


03 : 


59 


3 




2x2 


16 


Seeing 0.8 - 1.2 . High humidity 


5th Mar 




uU 


23 : 


52 


3 




2x2 


2 


Clear, but variable seeing up to 2.0" 


1st Jan 2009 


uu 




03 : 


06 


1.8 (sV); 3.6 


(«') 


lxl 


5 


Data gaps due to hardware problem. 




















Seeing 0.8" with some cloud 


2nd Jan 




ten 


07 : 


09 


1.8 (sV); 3.6 


in') 


lxl 


16 


Initial poor seeing (1.0 - 2.0") improves to 0.8". 




















Fair transparency. 


3rd Jan 


no 

VIZ 


u i 


04 : 


50 


1.8 (sV); 3.6 (u') 


1 x 1 


3 


Seeing 0.8", Good transparency. 
















LT/RISE 




17 Feb 2009 


23 


08 


00 : 


08 


30 




2x2 


2 


Seeing 0.5", high humidity 


14 Mar 


on 


OO 


21 : 


53 


30 




2x2 


2 


Seeing 2 — 3" , photometric 


21 Mar 


23 


09 


00 : 


09 


30 




2x2 


2 


Seeing 2 — 3" , photometric 


30 Mar 


nn 
UU 


U f 


02 : 


07 


30 




2x2 


4 


Seeing 2", photometric 


31 Mar 


00 


56 


01 : 


56 


30 




2x2 


2 


Seeing 0.5", photometric 


12 Apr 


21 


09 


22 : 


09 


30 




2x2 


2 


Seeing 0.5", some cloud 


19 Apr 


21 


52 


22 : 


52 


30 




2x2 


2 


Seeing 2 — 2.5", photometric 


13 May 


21 


25 


22 : 


25 


30 




2x2 


2 


Seeing 2", photometric 
















XMM-Newton 




23 Nov 2006 


10 


53 


20 : 


50 


35.7 x 10 3 






21 


EPIC MOS & RGS detectors 




11 


15 


20 : 


50 


34.0 x 10 3 






20 


EPIC pn detector 




10 


57 


15 : 


37 


10 






10 


Optical Monitor, fast mode with UVW1 filter 




15 


37 


20 : 


51 


10 






11 


Optical Monitor, fast mode with UVM2 filter 



with a 2 x 2 binning and the RISE V + R filter. Each ob- 
servation was lh in length (except the observation on 30 
March, which was twice as long) with exposure times of 
30s. The purpose of these observations was to characterise 
the superhump, since this exposure time is too long to ad- 
equately sample the eclipse. We reduced these data using 
aperture photometry as with the WHT/ULTRACAM data, 
using the ULTRACAM pipeline software. 

2.3 XMM-NEWTON 

SDSS 0926 was observed with XMM-Newton on 23 Novem- 
ber 2006. It was observed for 34.0 ksec in the EPIC pn de- 
tector and 35.7 ksec in the EPIC MOS detectors. It was 
detected with a mean count rate of 0.033±0.001 ct/s in the 
EPIC pn and 0.021±0.008 ct/s in the EPIC MOS (1+2), but 
was too faint in X-rays to be detected in the RGS detectors. 
The particle/X-ray background was low during the course of 
the observation. The Optical Monitor (OM) was configured 
in fast-mode and the observation time was split between the 
UVW1 (2450-3200A) and UVM2 (2050-2450A) filters. The 
source was detected with a mean count rate of 0.68 ct/s in 
UVW1 and 0.24 ct/s in UVM2. 

The X-ray data were processed using the XMM-JVewton 
Science Analysis Software (SAS) v9.0. For the EPIC detec- 
tors, data were extracted using an aperture of 30 centered 
on the source position. Background data were extracted 
from a source free region. The background data were scaled 
and subtracted from the source data. The OM data were 
reduced using omf chain. 



2.4 The Catalina Real-Time Transient Survey 

The Catalina Sky Survey (jLarson e t al. 1998|) is a search for 
near-Earth objects using the 0.7m f/1.9 Catalina Schmidt 
Telescope north of Tucson, Arizona. This survey uses a sin- 
gle unfiltered 4k x 4k CCD with 2.5" pixels, giving an 8 
deg 2 field of view. The Catali na Real-Time Transient Sur- 
vey (CRTS; iDrake et al.ll2009h began analysing these data 
in real-time in November 2007 for optical transients. 

The CRTS dataset contains 202 separate observations 
of SDSS 0926 between 10 November 2004 and 11 June 2010. 
Each observation is 30s in length, and they are divided up 
into groups of (typically) 4 observations taken over a ^30 
min p eriod. These data were reduced following IDrake et al.l 
l|2009h . 



3 RESULTS 

3.1 WHT/ULTRACAM light curves 

The March 2006 data are plotted in Figure [T] and the Jan- 
uary 2009 data are plotted in Figure [2] Additionally, we 
phase-folded the data on a night-by-night basis using the 
ephemeris given in Section [4. II and plot the results in Figure 
[3] We omit from this plot the short section of data collected 
on 5th March 2006. 

If we first examine the 2006 data, it is apparent that 
there are gross differences in the light curve from night to 
night. The shape of the eclipse features remains the same, 
but the superhump precesses through the light curve, and 
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Figure 1. Light curves of SDSS0926, observed in March 2006 with WHT/ULTRACAM. All data were collected simultaneously in the 
u'- (top, blue), g'- (middle, green) and r'-bands (bottom, red). For clarity we apply offsets of 0.05mjy to the g'-band data and O.lmjy 
to the n'-band data. The gaps in the first night of data are due to poor weather conditions. 



so we see the peak of the superhump emission at different 
phases on different nights. On 1 March the peak of the su- 
perhump is soon after the eclipse. On 2 March it is shortly 
before the eclipse, and on 3 March it is not immediately 
apparent, but the shape of the light curve before and af- 
ter the eclipse suggests that the superhump and the eclipse 
are approximately superimposed. If we examine the eclipse 
feature itself, we see that the primary eclipse is immedi- 
ately followed by a distinct second, smaller eclipse (this is 
most apparent in the 3 March data). We will show in Sec- 
tion 14.21 that these two eclipses are of the white dwarf and 
the bright spot, respectively. The eclipses are preceded by a 
small orbital 'hump' caused by the bright spot moving into 
the field of view. This is not immediately apparent since the 
bright spot is relatively weak in these data, so the out-of- 
eclipse variation is dominated by the superhump. As well 
as the superhump and eclipse features we see the stochas- 
tic 'flickering' variation that is characteristic of accreting 



sources. This variation is mitigated to some d egree in the 
phase -folded lightcurves (Figure [3]). Following ISmith et al.l 
(2002) we find the mean magnitudes outside of eclipse to be 
19.05 ± 0.10, 19.24 ± 0.07 and 19.39 ± 0.08 in it', g' and r' 
respectively. 

In contrast to 2006, in the 2009 data the shape of the 
out-of-eclipse light curve is roughly constant from night to 
night: we do not see the large variations caused by a super- 
hump component precessing through the light curve. The 
shape of the light curve on all three nights is most similar to 
the 2 March 2006 data, with the peak of the emission shortly 
before the eclipse. Since the position of this peak does not 
vary from night-to-night it is most likely due to the bright 
spot, and thus there seems to be no significant superhump 
contribution in these data. The mean magnitudes outside of 
eclipse are 18.94 ± 0.13, 19.31 ± 0.07 and 19.43 ± 0.11 in u' , 
g' and r' respectively, consistent with the 2006 values. Note 
also that the out-of-eclipse variation is double-humped, with 



© 0000 RAS, MNRAS 000, 000-000 



SDSS J0926+3624-' The shortest period eclipsing binary star 5 



i ■ 1 — 

1st January 2009 




3rd January 2009 



0.2 
0.15 

0.1 
0.05 




4 5 
Time (h) 



Figure 2. Light curves of SDSS0926, observed in January 2009 with WHT/ULTRACAM. All data were collected simultaneously in the 
u'- (top, blue), g'- (middle, green) and r'-bands (bottom, red). For clarity we apply offsets of 0.05mjy to the (/-band data and O.lmjy 
to the n'-band data. The gaps in the first night of data are due to poor weather conditions and a hardware fault. 



a peak at a phase of ~0.3 as well as the main peak at ~0.8. 
The most likely explanation for this is that bright spot is ver- 
tically extended, or disc is optically thin, so we are seeing 
emission from the bright spot all the way round the orbit. 



3.2 Non-orbital variability in the ULTRACAM 
data 

In Section [4] we use the high time resolution ULTRACAM 
data to make precise parameter determinations. However, 
first it is necessary to examine the non-orbital variabil- 
ity in this system. We begin by examining quasi-periodic 
variability in the ULTRACAM data. Secondly, it is impor- 
tant to characterise the superhumps present in the 2006 
observations, so these features can be subtracted from the 
lightcurves. 

3.2.1 Quasi- coherent variability 

In Figure[4]we plot Lomb-Scargle periodograms (|Press!l2002l ) 
for the complete 2006 and 2009 (/-band datasets obtained 
with WHT/ULTRACAM. In the 2006 data we detect a sig- 
nal at a frequency of ~1700cycles/day (P~50s), although it 
is incoherent and spread over a wide frequency range. We 
estimate the quality factor (the peak centroid frequency di- 
vided by its full width at half maximum) of this signal to be 
Q ~4 in the (/-band data. We computed the periodograms 



for each of the first three 2006 nights separately and we 
detected this signal every night. The signal is high in the g'- 
band, and is barely detected in the v! band. There is possibly 
a signal at the lower frequency of ~1400 cycles/day in the 
2009 data, but it is much weaker than the 2006 signal. We 
did not find any signals in the periodograms at higher fre- 
quencies beyond the 5000 cycles/day range plotted in Figure 

a 

Similar quasi-coherent var iability was first observed 
in CVs some deca des ago (jWarner fc Robinson! Il972l ; 
iPatterson et al.ll 19771 ) , and ha s since been observed in many 
CVs and X-ray binaries (see IWarner fe Woudtl [2008 for a 
recent review). The peak period of the signal we detect in 
the 2006 data is low for a QPO, but this may be due to the 
fact that the accretion disc in an ultracompact binary such 
as SDSS 0926 is much smaller and less massive than the disc 
in conventional CV systems. 



3.2.2 Superhumps 

Figure [5] shows (/-band Lomb-Scargle periodgrams in the 
vicinity of the orbital frequency, with various manipula- 
tions applied. The top panel shows the combined 2006 
dataset, with no modification. The second panel uses the 
same dataset, but with the eclipse features masked. A phase 
range of 0.18 is masked, centred on the mid-point of the 
white dwarf eclipse. This phase range is sufficient to cover 
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Figure 3. Phase folded and binned light curves, showing the superhump variation from night-to-night. In the top row we plot the first 
three nights of data collected in 2006. In the bottom row we plot the three nights of data taken in 2009. We plot separately the data in 
the u'- (top, blue), g'- (middle, green) and r'-bands (bottom, red). 
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Figure 4. Lomb-Scarglc periodgrams for the g'-band 2006 (top) 
and 2009 (bottom) WHT /ULTRAC AM datasets. We convert 
both the y- and x-axes to a logarithm scale, and then uniformly 
bin the data along the x-axis (frequency). In the 2006 dataset we 
see a quasi-periodic oscillation (QPO) at a frequency of ~1700 
cycles/day (P=~50s). 



Table 2. Superhump parameters. We list the period and am- 
plitude of the primary sine function fitted to the complete 2006 
dataset, for each of the three bands. We list also the period excess 
e = (Pg n — Porb) I Porb, where Psh an d Porb are th e superhump 
and orbital periods, respectively. 

Filter Period Amplitude Period excess (e) 

(min) (mJy) 

u' 28.560 ±0.003 44.1 ±0.8 0.00875 ± 0.00010 

g' 28.558 ±0.001 50.6 ±0.5 0.00869 ± 0.00005 

r' 28.553 ±0.001 60.1 ±0.5 0.00849 ± 0.00005 



the eclipses of both the white dwarf and bright spot and the 
peak of the bright spot emission. The third panel shows the 
same dataset with the eclipse features unmasked but with 
the superhump subtracted (see below). The bottom panel 
shows the combined 2009 dataset, with no modification. In 
all four panels, vertical lines mark the positions of the su- 
perhump and orbital frequencies. 

If the top panel of Figure [5] is examined, it can be 
seen that the power due to superhumps is clearly apparent, 
stronger than the orbital signal and peaking at a slightly 
lower frequency. The two signals are confused in this first 
panel, but the superhump is seen as being clearly distinct 
in the second panel, in which the majority of the orbital 
modulation is masked. 

In order to determine the parameters of the 2006 super- 
hump, we fitted a model to the combined dataset with the 
eclipse features masked. The model consists of a combination 
of six 4-parameter sine functions: three for the superhump, 
fitting the primary frequency and the second and third har- 
monics, and three for harmonics of the orbital period, so as 
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Figure 5. Lomb-Scargle periodograms for the g' band data. We 
plot the power around the orbital frequency, with the peak of each 
distribution normalised. The vertical lines mark the frequency of 
the superhump (blue, dotted) and the orbital (red, solid) mod- 
ulations. The four periodograms are (from top to bottom): (i) 
The 2006 data, showing both the superhump and orbital modu- 
lations, (ii) The 2006 data with the eclipse and bright spot fea- 
tures masked, (iii) The 2009 data with the superhump subtracted, 
using the method described in Section 13.2.21 (iv) The 2009 data. 



to fit any residual signal left after the masking of the eclipse 
features. In the third panel of Figure [5] the periodogram for 
the unmasked dataset with the superhump components sub- 
tracted is plotted. We see that our model fits do a good job 
of cleaning the superhump signal from the data. 

Our findings for the superhump period are listed in 
Table [2] The uncertainties on these periods were deter- 
mined from fits to a series of sample data sets derived 
from the originals using the bootstrap method |Efronlll979l ; 
lEfron fc Tibshiranilll993l ). The amplitudes of the superhump 
harmonics are < 10% of the amplitude values for the pri- 
mary frequency. We find the period of the variation to be 
consistent at the lcr level for the u'- and p'-bands, but the 
r'-band period is lower. This inconsistency is probably due 
to some extra intrinsic variability, such as accretion-driven 
flickering. The amplitude of the modulation increases at 
longer wavelengths. We list also in this table the period 
excesses e in each band, using the orbital period given in 
Section I4T1 

Finally, the fourth panel of Figure [5] shows the 2009 
data. This plot shows a clean signal at the orbital frequency, 
confirming that there is no superhump modulation in these 
data. Following these observations, we obtained a series of 
lh light curves using LT/RISE over the first half of 2009. 
The purpose of these observations was to examine the long 
term behaviour of the superhump, although with hindsight 
they were perhaps too short. Two of these light curves were 
in the immediate aftermath of an outburst in this system, 
and we will discuss these separately in Section 13.3.11 Of the 
remaining light curves, few show clear evidence for the su- 
perhump. On 17 February and 21 March there is a 'hump' 
just before the eclipse, but this could be due to the bright 
spot. Most of the remaining quiescent light curves show lit- 
tle out-of-eclipse variation. One exception is the light curve 
obtained on 19 April, which shows a clear brightening imme- 
diately following the eclipse, which can only be explained by 
the superhump. This light curve was obtained 21 days after 
the detection of the outburst in this system. There was no 
clear evidence for the superhump obtained seven days prior 
to this one, on 12 April. 

IPatterson et all (|2005h suggested e = 0.18q + 0.29g 2 as 
an empirical relationship between the superhump period ex- 
cess e = (Psh — Porb) / Port, and the mass ratio q. This 
relationship was calibrated using measure ments of a series 
of ecl ipsing systems, listed in Table 7 of IPatterson et al.l 
(2005). The relation is pinned by setting e = when q = 0, 
but this is an assumption and not empirically determined: 
of the calibration systems only KV UMa has a mass ra- 
tio < 0.05 and this determination is very uncertain, and 
so the calibration is potentially poor at low mass ratios. 
SDSS 0926 is therefore potentially a strong test of this re- 
lationship, although it has been argued t hat it would not 
apply to AM CVn syste ms ( Pearson! 2007 ) . In Figure [6] we 
reproduce Figure 1 from IPatterson et al T |2005l ). adding our 
measurements from the 2006 data in u , g' and r' , using 
the values given in Tables [2] and [6] As well as the Pat- 
terson relati on, we also plot th e slightly modified relation 
proposed by iKato et al.l (|2009l ) (e = 0.16(2 ) <? + 25(7)q 2 ) 
and the linear relation proposed by iKnigge! (|2006h (?(e) = 
(0.114±0.005) + (3.97±0.41) x (e-0.025)). The Knigge rela- 
tion does not assume e = when q = 0. Our measurements 
of the period excess in SDSS 0926 are consistent with all of 
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Figure 6. Mass ratio q versus the superhump period excess e. 
The blue, green and red points show our u'-, g'- and r'-band 
determinations, respectively. The region around these points is 
magnified in the inset. The solid line is the e = 0.18q + 0.29g 2 
relationship proposed bv lPatterson et"aL I te005|V The dashed line 
is the e = 0.16<j + 0.25<j 2 relation proposed bv lKato etahl J2009T) . 
The dotted line is the linear q(e) = (0*114±0.00 5) + (3.97±0.41) X 
(e — 0.025) relation proposed by Kniggc (200j]). The black points 
are the eclipsing CVs listed as calibration sources in Table 7 of 
IPatterson et alj (2005h . 

these r elations to within their uncertainties, and are closest 
to the IPatterson et alj (|2005l ) relation, which suggests the 
assumption of e = when q — is reasonable. 

3.3 Outbursting behaviour 

The superhumps are a transient phenomenon that is driven 
by outbursts, so here we examine the outbursting history of 
this system, presenting the first observations of SDSS 0926 
in outburst. We begin by discussing observations taken with 
the LT in the immediate aftermath of an outburst in March 
2009. We go on to examine the long-term outbursting be- 
haviour using five years of CRTS observations. 

3.3.1 The March 2009 outburst 

On 29 March 2009 it was discovered that SDSS 0926 was 
in outburst. We obtained 2h of data with LT/RISE on the 
subsequent night, and a further lh on the night after that. 
We plot these data in Figure [7] The flux scale for these data 
is such that the mean, quiescent out-of eclipse flux is equal 
to 1. We also plot the model fit (described in Section 14.21 
below) to our quiescent superhump subtracted LT/RISE ob- 
servations for comparison. We apply an arbitrary flux offset 
to this model light curve so as to overlay it on the outburst 
data. We see that at the beginning of our 30 March obser- 
vation the system was a factor of ~3.5 brighter than during 
quiescence and declining quickly: this has dropped to a fac- 
tor of 3 by the end of the 2h observation. The following night 
this has dropped to a factor of 2. There is still some evidence 
of the outburst in our next observation on 12 April: the out- 
of-eclipse flux in this light curve is ~10% greater than the 
mean level. 

We phase- folded the 30 March data using the ephemeris 
given in Section [4.11 and plot the results in Figure [8] Again, 



we plot a quiescent model light curve for comparison. In 
this plot the structure in the light curve after the outburst 
is more evident. There are two main features: the primary 
eclipse centred on a phase of 0, and a separate, smaller dim- 
ming in the light curve centred on a phase of ~— 0.25. If we 
examine the primary eclipse first we see it is much deeper 
and wider than the quiescent white dwarf eclipse. The eclipse 
width is consistent with this being an eclipse of the accre- 
tion disc. The eclipse is asymmetric however, and there is a 
clear 'step' in the egress. This suggests an uneven distribu- 
tion of flux over the surface of the accretion disc itself, or 
(perhaps more likely) this step could be due to the bright 
spot. The step height of this egress is much larger than the 
bright spot egress during quiescence however, and so if this 
feature is attributed to the spot it would imply an enhanced 
mass transfer during the outbursting state, or an increased 
viscous heating at the bright spot position. 

The second unusual feature in the phase-folded light 
curve is the dimming centred on a phase of ~— 0.25. This 
dimming begins at approximately the same time that the 
bright spot begins to come into view in the quiescent light 
curve. The phase of this feature is such that the eclipsing 
component cannot be either the white dwarf, the donor or 
the bright spot. We su ggest this fea ture may be due to a 
warped accretion disc jpringld[l996h , the distortion being 
induced by the outburst. The dimming we see can therefore 
be explained by obscuration of the white dwarf or inner disc 
by the accretion disc itself. This is a very short lived fea- 
ture: we see no evidence for it in the data collected on the 
subsequent night. 

Finally, as we discussed in Section 13.2.21 we detect su- 
perhumps in this system in data obtained 21 days after the 
outburst. This is consistent with the expectation that, since 
the superhumps in this system are not permanently present, 
they are induced by the perturbation of the disc by the out- 
burst. 

3.3.2 Long-term variations 

We plot in Figure [9] the unfiltered data collected by the 
Catalina Real-Time Transient Survey (CRTS). These data 
were collected between 10 November 2004 and 11 June 2010. 
We observe a total of either six or seven outbursts in these 
data, but the time sampling is such that they do not give an 
exhaustive picture of the outbursting behaviour of the source 
over this time period: there are no observations around the 
time of the March 2009 outburst, for example, and so this 
outburst is missed. Note that we do detect an outburst ~25 
days before the 2006 ULTRACAM observations. It is likely 
that this outburst induced the superhumps we see in these 
data by perturbing the disc. 

The CRTS data is split up into blocks of four 30s ob- 
servations taken over a ~30 min period. There is typically 
> 20 days between these blocks, and so in most cases the 
source has returned to quiescence by the time of the observa- 
tion block subsequent to the detection of the outburst. The 
exception is the sixth and seventh outburst: these are sepa- 
rated by only eight days with no intermediate points, so may 
be one long outburst or two separate ones. An 8+ day out- 
burst is long considering the orbital period of this system: 
if we scale by orbital period this would be the equivalent of 
a 2 month long outburst in a P = 3.5 h dwarf nova. How- 
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Figure 7. LT/RISE light curves obtained on 30 and 31 March 2009: one day and two days after an outburst in this source was first 
detected, respectively. We plot these data in flux units, with the mean quiescent out-of-eclipse flux set to unity. The solid red line is a 
model fit to our quiescent, superhump-subtracted LT/RISE observations with an arbitrary flux offset, which we plot for comparative 
purposes. 
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Figure 8. The LT/RISE data collected one day after the March 
2009 outburst, phase-folded using the ephemeris given in Section 
I4.ll We plot these data in flux units, with the mean quiescent 
out-of-eclipse flux set to unity. The solid red line is the quiescent, 
superhump-subtracted model light curve with an arbitrary flux 
offset, which we plot for comparative purposes. 



ever, outbursts of ~12 days hav e been observed in K L Dra, 
a P = 25 min AM CVn system ijRamsav et alJl201Cl ). 

The mean magnitude of the quiescent, out-of-eclipse 
points is 19.33 ± 0.31. The outbursts range in magnitude 
from 17.11 ± 0.08 to 16.81 ± 0.12. Note that for all six out- 
bursts it is unclear how much time has elapsed between the 
initial rise and the observation, and so these measurements 
set a lower limit on the peak brightness. Our LT/RISE ob- 
servations of the March 2009 outburst showed a factor ~3 
enhancement in flux a day after the initial detection, and a 
factor ~2 enhancement a day after that. This is relatively 
modest compared to the outbursts observed in the CRTS 
data. We conclude that the March 2009 outburst was a rel- 
atively minor one for this source. An outburst of this size 
would only be detectable in the CRTS data within 1-2 days 
of the initial rise. The average outburst recurrence time is 
therefore difficult to determine. If we exclude the two out- 
bursts which occur within 8 days of each other the time be- 
tween observed outbursts ranges from 104 to 449 days, but 
since we may have missed a number of outbursts the actual 
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Figure 9. Observations made by the Catalina Real-Time Tran- 
sient Survey (CRTS). These data were collected between 10 
November 2004 and 11 June 2010, and six outbursts were ob- 
served over that period. The solid lines at the top of the plot 
indicate the times of the 2006 and 2009 observations with the 
WHT and LT. 



recurrence time may be shorter. For comparis on, an outburst 
cycle of ~60 days was observed in KL Dra ijRamsav et al.l 
|2010| ). 



3.4 XMM-Newton observations 

Observations have sho wn that the ultravio let luminosities 
of AM CVns are high (|Ramsav et alJ [ 20051 ). hence to fully 
understand the energy budget of these systems it is neces- 
sary for us to make observations at X-ray and ultraviolet 
wavelengths. We introduce here the data we obtained with 
XMM-Newton, and go on to use the UV OM data to deter- 
mine the temperature of the primary white dwarf. 
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N H 



k'1 77 



1.3±0.3 x 10 21 cm" 2 

1 - 8 -0.5 

30 keV (fix) 

h2.8 w in-13 



Observed flux 0.1-10keV 1-21q i X iO -13 er g s_1 cm ~ 2 
Unabsorbed 

2.7±Q g X 10~ 13 erg s^ 1 cm- 2 



bolometric flux 
Unabsorbed 
bolometric luminosity 



3.2+8-j x 10 29 erg s' 1 d 2 
X 2 =l-03 (103 dof) 



Table 3. The X-ray spectral fit parameters derived using a si- 
multaneous fit to the EPIC pn and MOS spectra. 



3.4-1 X-ray data 

We extracted light curves from the EPIC pn and both EPIC 
MOS cameras and combined them into one light curve, 
which was folded using the ephemeris given in Section 14.11 
We found no evidence for an eclipse in the X-ray light curve. 
To investigate this further we searched for periods using a 
Discrete Fourier Transform and phase dispersion methods 
and found no clear signal of an eclipse. 

Using the original X-ray light curve as a benchmark, we 
added an eclipse of given depth to the light curve. We gen- 
erated 100 light curves for a given eclipse depth. We found 
that for a partial eclipse with a depth of < 90%, we would 
likely not detect the eclipse. Even for a total eclipse we would 
have only a 70% probability of detecting the eclipse. 

We extracted integrated X-ray spectra from the EPIC 
pn and MOS detectors (with corresponding background 
spectra) and fitted them simultaneously using an absorbed 
thermal plasma model. Unlike other AM CVn systems 
ijRamsav et al.ll2005l . 2006), we found that varying the metal 
abundance did not significantly improve the goodness of fit 
(probably since the signal-to-noise ratio of the spectra was 
relatively low). We obtained a good fit to the data, and the 
derived spectral parameters are listed in Table [3] 

The best-fit value for the absorption is high (1.3±0.3 x 
10 21 cm -2 ) for an object at a high galactic latitude (+46°). 
In contrast, the abs orption to the edge of th e Galaxy is 
~ 1.4 x 10 20 cm~ 2 l|Dickev fc Lockmanl I1990D . If we fix 
the column density para meter in the spectral fits to the 
iDickev &: Lockmanl (|l990l ) value the resulting fit (x^=l-96, 
104 dof) was poorer at a confidence level of > 99.99%. If 
this 'extra' absorption originates in SDSS 0926 itself it may 
be due to viewing the boundary layer through the disc, since 
the binary inclination is high. Finally, we calculate the unab- 
sorbed bolometric X-ray luminosity using the distance given 

3 —1 
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Figure 10. Phase folded UV lightcurves obtained using XMM- 
Newton OM. The top and bottom plots are through the UVW1 
and UVM2 filters, respectively. 



3-4-2 Ultraviolet data 

We binned the OM light curves into 5s bins and folded them 
on the ephemeris given in Section 14.11 The resulting light 
curve is plotted in Figure fTCTl The eclipse is clearly seen in 
the UVW1 filter, where the count rate is consistent with 
zero at 0=0.0 in a single bin. In the shorter wavelength fil- 
ter, UVM2, the signal-to-noise ratio is lower, but there is 
some evidence for a dip at 0=0.0, although it is not to- 
tal. We searched for other periods using a Discrete Fourier 
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in Section 13.4.11 the galactic E(B — V) is negligible, but 
we measure some extra absorption in the X-ray data. This 
is probably related to the high inclination of this system 
and due to our looking through a significant amount of gas 
above the accretion disc. Si milar effects are seen in a number 
of CVs, such as V893 Sco (jMukai et alJl2009l ). This mate- 
rial will not be in the form of dust and so will not cause 
'extinction' in the classical sense, but it is possible that it 
may caus e an 'accretion cu rtain' in this system, such as in 
OY Car l|Horne et alJll994h . If this was the case then our 
temperature determination was a lower limit. However, it is 
impossible to detect such an effect without UV spectroscopic 
data. Our temperature determination for the white dwarf 
is consiste nt (within the uncert ainties) with the theoretical 
models of Bi ldsten et all (|2006l ). which predict a tempera- 
ture of 18, 000K for this system. 



Figure 11. We determine the white dwarf temperature by fit- 
ting our measurements of the white d warf flux in the optic al and 
UV to a library of synthetic spectra (Gansickc et al. 1995). The 
spectrum we plot here is for a 17, 000K white dwarf, and is the 
best fit to our flux determinations. The numbers in the top right 
corner are the temperature and log g for the model fit. 



Transform toeemindll975h and phase dispersion methods 
|Stellingweri 19781 ) and found no other significant periods 
in the light curves. 

We derived UV fluxes by converting the count rate in 
the two UV filters to a flux by assuming a conversion fac- 
tor which was derived using observations of white dwarfs 
using XMM-Newton OM. A mean count rate of 0.71 ± 0.02 
ct/s in the UVW1 filter corresponds to a flux of 3.2xl0~ 16 
erg s _1 cm" 2 /A, while 0.27 ± 0.01 ct/s in the UVM2 filter 
corresponds to a flux of 5.3x 10~ 16 erg s _1 cm~ 2 /A. We cal- 
culate the ult r aviole t lumi nosity following the approach of 
iRamsav et al.l l|2005l . l2006t ). by assuming a blackbody flux 
distribution and fixing the normalisation to give the inferred 
de-reddened ultraviolet flux in the two filters. This gives an 
ultraviolet luminosity of L uv = 1 — 3 X 10 32 erg s _1 . Both 
this luminosity and the X-ray luminosity L x , given in the 
previous section, are consistent wit h the findings for ot her 
AM CVn systems of similar period ()Ramsav et al.ll2006h . 

We used the UV and optical fluxes to determine the 
effective temperature of the white dwarf, by estimating the 
contribution of the white dwarf in each band. In the three 
ULTRACAM bands this parameter is determined through 
our model fits, as discussed in Section [4] In the UV data 
we estimated the white dwarf contribution from the eclipse 
depth. We made a correction to compensate for the fact 
that the eclipse is only partial, by determining the frac- 
tion of the white dwarf's surface area that is obscured dur- 
ing the eclipse in our optical fits. This correction is small 
though, as evidenced by the fact that the UVW1 flux reaches 
zero during the eclipse, within the errors. We calculated the 
white dwarf contributions to be UVM2 ~0.04±0.02, UVW1 
~0.062 ± 0.012, v! = 0.057 ± 0.005, g' = 0.046 ± 0.002 
and r' = 0.035 ± 0.002 mjy. We then fitted these fluxes 
with the whit e dwa rf model atmospheres introduced in 
iGimsicke et ail (|l995l ). We find a white dwarf temperature 
of 17, 000K to be consistent with our measurements (Fig- 
ure [11]). Note that in making this estimate we have not ap- 
plied an extinction correction to these fluxes. As we noted 



4 LIGHT CURVE ANALYSIS 

In this section we describe the model we fitted to the phase- 
folded and binned WHT/ULTRACAM 2006 and 2009 data 
in order to make parameter determinations for this system. 
For the 2006 data it was first necessary to subtract the night- 
to-night variations caused by the superhump; this is detailed 
in Section [3. 2. 2 1 Secondly, we fitted the data on an eclipse- 
by-eclipse basis in order to refine the ephemeris, which we 
discuss in Section [4. II We then phase-folded the data using 
this ephemeris and fitted it with a Markov Chain Monte 
Carlo (MCMC) light-curve model in order to obtain final 
parameter determinations. We fit the three bands separately. 
This fitting process is described in Section [4. 21 



4.1 Determination of the ephemeris 

We divided each light curve into separate orbital cycles, and 
then fit the model determined in Section [4. 2 | to each cy cle 
using the Levendburg-Marquardt method ( Press! I2002T I in 
order to generate a series of eclipse timings. These timings 
are on the barycentric dynamical timescale, corrected for 
light travel to the solar system barycentre. A least-squares 
fit to all of these data yields the ephemeris 



MJD(TDB) = 53795.9455191(5) + 0.01966127289(2)^ 

for the mid-point of the white dwarf eclipse. We plot the 
residuals of the linear ephemeris in Figure [121 There is some 
systematic variation, perhaps due to the residuals of the su- 
perhump or flickering. With only two epochs of observation 
any long-term departure from a linear ephemeris cannot be 
determined: a third epoch of high-speed observation will be 
necessary to identify any period changes in this system. 



4.2 Fitting the phase-folded light curves 

In order to make precise parameter determinations we chose 
to combine our data into phase-folded and binned light 
curves. We began by preparing the 2006 data for fitting by 
subtracting the superhump modulation from the data, us- 
ing the parameters determined in Section 13.2.21 Since the 
superhump is not seen in 2009, this step is not necessary for 



© 0000 RAS, MNRAS 000, 000-000 



12 CM. Copper-wheat et al. 




2 4 6 

3 
2 

_ 1 
o o 

° -1 

-2 
-3 



10 12 14 46 48 50 52 54 56 58 60 62 64 



96 98 100 102 104 106 108 110 112 114 204 206 




52744 52746 52748 52750 52752 52794 52796 52798 52800 52802 52804 52806 52808 52810 52812 52854 52856 52858 

Cycle 

Figure 12. (O — C) values plotted against cycle number, using the eclipse timings determined from the model fits to our 
WHT/ULTRACAM data and the linear ephemeris given in Section l4.f I Cycle corresponds to the first eclipse in the 1 March 2006 
dataset. Note the break in the x-axis of this plot between the various nights of data. 



this second epoch of data. Once the superhump was sub- 
tracted it became clear that other binary parameters var- 
ied over the course of our observations, in particular the 
disc radius which changes significantly from night to night. 
These variations prevented us from combining our entire 
WHT/ULTRACAM dataset. We therefore chose to create 
separate phase-folded light curves for each individual night 
of observation. Even this could potentially introduce some 
systematic effect on our results due to the change in disc 
radius between the beginning and end of each night's obser- 
vation. We examine this in detail in Section [5~4l and find that 
such effects are small, so we do not believe this influences 
our results to a significant degree. Since we cannot com- 
bine nights, we chose to omit the nights of 5 March 2006, 
1 January 2009 and 3 January 2009, in which we only have 
a few cycles of data. We therefore fitted a total of twelve 
light curves: the u'-, g'- and r'-band phase- folded data for 
the nights of 1, 2 and 3 of March 2006; and 2 January 2009. 

We modelled the light curve with LCURVE, a code de- 
veloped to fit light-curves characteristic of eclipsing dwarf 
novae and detached white dwarf / M dwarf binary stars. A 
co mplete description of this code is given in the appendix 
of ICopperwheat et al.l (|2010t ). We implemented this code in 
this work with two modifications, both to the bright spot 
component. In lCopperwheat et al.l (|2010t ) the bright spot is 
modelled as a line of elements which lie upon a straight line 
in the orbital plane. The surface brightness of the elements 
is parameterised with two power-law exponents. Since the 
bright spot in SDSS 0926 is a relatively weak component 



of the emission we do not require this degree of complexity. 
We therefore use a simpler bright spot model, setting the 
exponent 7 to 1 (our bright sp ot model is therefo re identical 
to the earlier prescription o f iHorne et al1 ll994l). Secondly, 
in lCopperwheat et al.l (|2010T ) the angle 4> was defined as the 
angle the line of elements of the bright spot makes with the 
line of centres between the two stars. We have changed the 
definition of this angle in our code: <j> is now the angle the 
line of elements of the bright spot makes with the tangent 
to the outer edge of the accretion disc. This modification 
makes this angle easier to interpret in a physical context, 
since 4> = implies a bright spot which runs along the outer 
edge of the disc, and so we would expect </> to be close to 
this. We noted in ICopperwheat et al.l l|2010l ) that the two 
angles (f> and ip (the angle away from the perpendicular at 
which the light from the bright spot is beamed) are highly 
correlated and tend to be poorly constrained, and so here 
we fix <f> — 0. 

In the LCURVE code the binary is defined by four com- 
ponents: a white dwarf primary, a Roche-lobe filling sec- 
ondary star, accretion disc and bright-spot. We first ob- 
tained an initial fit to each light curve using the simplex and 
Levenberg-Marquardt methods l|Pressl 2002). We then used 
a Markov Chain Monte Carlo (MCMC) algorithm for min- 
imisation and determination of u ncertainties (details of ou r 
MCMC method are also given in lCopperwheat et al.l2Qlol). 
We fit ted all of the parameters used in ICopperwheat et all 
(2010), with the exception of 7 and <j>, as described above. 
We additionally set the accretion disc radius Rdisc equal 
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Table 4. Limb darkening coeffic ients fo r the p rimary white dwarf, 
using the four-coefficient law of Claret ( 200(j) • 



Coefficient u' g 1 r' 



1 


1.35 


1.18 


1.37 


2 


-1.30 


-1.13 


-2.07 


3 


0.90 


0.75 


1.96 


1 


-0.28 


-0.23 


-0.72 



to Rspot'- the distance between the bright spot and the pri- 
mary white dwarf. We therefore assume that the 'head' of 
the bright spot is on the outer rim of the disc. Rdisc is 
ra ther poorly constra i ned b y our data, and as we showed 
in lCopperwheat et al.l l|2010l ) R sp ot is highly correlated with 
Rdisc, so this is a reasonable approximation. One additional 
parameter which our results are potentially sensitive to is 
the limb darkening of the primary white dwarf. We made 
an initial fit to the data in order to estimate the effective 
temperature and surface gravity (we discuss these quanti- 
ties in Sectio n | 5.2l). and the n used a model atmosphere code 
(described in lKoesterll2010l ) to calculate the specific inten- 
sity at different points across the stellar disc. We then fitted 
these values to determine the limb darkening coefficients. 
We tried various limb darkening la ws, and f ound the best fit 
was for the four-parameter law of IClared (|200Ch . although 
the choice between this and a fourth order polynomial is 
unlikely to influence our results. We list our determinations 
of the coefficients in Table [4] These values were used in all 
our MCMC fits. 

The parameter determinations from these fits are listed 
in Table [S] Note that the uncertainties on these MCMC re- 
sults are non-Gaussian, and so the values we quote in this 
table only provide an approximate description of the un- 
certainties. We plot the distribution of the mass ratio and 
inclination values from our fits in Figure 1141 Additionally, 
in Figure[l3]the phase-folded light curves for the four nights 
are plotted, along with the best model fits. If we examine 
Figure [131 first, we see we obtain a good fit to the data in 
all three bands for each of the four nights. The dominant 
component in the light curve is the primary white dwarf, 
and the main eclipse in the light curve is the eclipse of this 
feature. Note that the eclipse is round-bottomed, meaning 
it is a partial eclipse of the white dwarf. The second, smaller 
eclipse in these light curves is of the bright spot, which is 
a much weaker component. The phase of this eclipse clearly 
varies with respect to the phase of the white dwarf eclipse: 
this can be understood in terms of a change in the relative 
position of the bright spot due to variations in the disc ra- 
dius. Compare for example the first night of 2006 with the 
second. On the first night, the two eclipses are clearly dis- 
tinct, with the egress of the white dwarf followed by the 
ingress of the bright spot. On the second night, the egress 
and the ingress overlap. This implies a larger apparent ac- 
cretion disc radius on the second night, due to the precession 
of the elliptical disc. The two other nights of data lie some- 
where between these two extremes. The fact that the two 
eclipses are sometimes separate implies a very low mass ra- 
tio for this system: in higher mass ratio systems we would 
expect to see the ingress of both the white dwarf and the 
bright spot before the white dwarf egress. 

We now turn to the parameter determinations listed in 



Table [5] The first four of these are the parameters that are 
important in characterising the system: these are the mass 
ratio q, the binary inclination i, the primary white dwarf 
radius scaled by the binary separation Ri/a, and the ac- 
cretion disc radius scaled by the binary separation Rdisc/ a,. 
The mass ratio and the inclination are highly correlated, 
and we discuss these parameters in Section [5.11 The white 
dwarf temperature is used with q, i and Ri/a to derive the 
remaining binary parameters in Section [5.21 We discuss the 
disc radius in more detail in Section f5. 41 The remaining six 
parameters in Table [S] pertain to the accretion disc and the 
bright spot. We find these parameters to be rather poorly 
constrained in general, due to the fact that the disc and the 
bright spot make a relatively weak contribution to the flux in 
this system. However, we find no correlation between these 
parameters and the 'important' parameters listed above, and 
so the uncertainty in these values does not imply any further 
systematic uncertainty in the determinations of Sections l5.ll 
and El 



5 DISCUSSION 

5.1 Mass ratio and inclination 

If we assume a Roche-lobe filling donor star, the phase width 
of the white dwarf eclipse is then an observable quantity 
that is intrinsically linked to two physical properties: the 
mass ratio and the binary inclination. For a higher binary 
inclination the duration of the eclipse will be greater, thus 
to maintain the same phase width, as the inclination is in- 
creased, the size of the donor, and hence the mass ratio, 
must be decreased. There is th erefore a uni que relationship 
between these two properties (|Bailevl lT979). This degener- 
acy can be broken since we have an additional geometric 
constraint due to the ingress and egress of the bright spot. 
The path of the accretion stream and hence the position of 
the bright spot is modified by the mass ratio. With this ad- 
ditional information we can determine both the mass ratio 
and inclination in this system. 

Our MCMC results for the mass ratio versus inclination 
are plotted in Figure 1141 We see that these points are dis- 
tributed along a curved path across this plot: this is the line 
of constant phase width for the white dwarf, which we mea- 
sure to be ~0.0220. The points are constrained to this line, 
with the scatter due to the uncertainty in the phase width. 
We calculate the weighted means of the q, i and R\ja values 
in Table [5j and present the results in Table [6] We emphasise 
at this point that the weighted means may underestimate 
the uncertainties in these parameters. The MCMC results 
in Table [5] show variations which are formally significant, 
most likely due to effects such as flickering and residuals 
from the superhumps. In addition, the bright spot compo- 
nent is weak in this system, and so the best fit model for this 
component can be quite different in different bands, partic- 
ularly in u . It is important therefore to note that the results 
we discuss in this and the following sections, and in Table [6j 
may not fully account for these systematic effects, and the 
magnitude of these effects are better understood by referring 
to the individual night's results in Table [5] 

We find the mass ratio Ms/ Mi to be q = 0.041±0.002 in 
our (/-band measurements. The inclination is i = 82.6 ± 0.3 
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Figure 13. Phase-folded and binned light curves for the first, second and third nights of the 2006 observations, and the second night 
of the 2009 observations. For each we plot the three bands separately (top, u'; middle, </; bottom, r'). We plot the average flux in mjy 
against the binary phase, where a phase of corresponds to the mid-eclipse of the white dwarf. We plot the datapoints with uncertainties 
in black, and the best model fits to these data in blue, green and red (for u' , g' and r', respectively). For clarity we apply offsets of 
0.005mjy to the g-band data and O.OlmJy to the u-band. For each night we also show underneath the lightcurve the three components 
of the g'-band model plotted separately, showing the relative strengths of the bright spot, accretion disc and white dwarf. These three 
lines have also been offset for clarity. 



deg in g' . The values of q and i in the other two bands are 
consistent w ith these values. Ou r 2006 data were originally 
published in iMarsh et all {2007) and there we reported dif- 
ferent values of 0.035 ± 0.002 and 83.1 ± 0.1 deg for q and 
i. This discrepancy is due to a number of factors. Primar- 
ily, the MCMC method we use for minimisation is supe- 
rior to the L e venbe rg-Marquardt minimisation employed in 
IMarsh et all l|2007l ). We have found in situations of strong 
degeneracy, such as between q and i here, the Levenberg- 
Marquardt and simplex methods have a tendency to sto p 
before finding the true minimum |Copperwheat et alfe oiO). 
Additionally, we have used more appropriate limb-darkening 
coefficients in this work, since we have been a ble to estimate 
the white dwarf temperature. IMarsh et alj (|2007l ) also ar- 
rived at their parameter determinations by combining data 
from different nights, which as we have discussed may intro- 



duce some systematic uncertainty due to the changing disc 
radius. 

5.2 Other binary parameters 

In Table [5] we list our determination of the primary white 
dwarf radius, scaled by the binary separation. We used this, 
combined with our determination of the mass ratio and the 
binary inclination, to calculate the remaining binary param- 
eters. One additional piece of information that was needed 
for this is a mass/radius relation for the primary white 
dwarf. 

We determined this by using the white dwarf tempera- 
ture, which we found in Section \3. 4. 21 to be 17, 000K. How- 
ever this is an iterative process, since the temperature de- 
termination itself requires a mass/radius relation. The pro- 
cess was as follows. We began by assuming a mass/radius 
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Figure 14. The distribution of the mass ratio q and inclination i results from our MCMC runs, with the green and red regions indicating 
the 1(7 and 3<r confidence intervals. We plot (from left to right) nights 1, 2 and 3 of the 2006 observations, and night 2 of the 2009 
observations. For each night we plot (top to bottom) the distribution in the u' , g' and r' bands. 



relation for the white dwarf, and for this we used the 
Eggleton zero-temperature m ass/radius relation (quoted in 
IVerbunt fc Rappaporj[l988l ). We subsequently determined 
the mass and radius of the white dwarf, and from our 
model calculated the white dwarf contribution in each band. 
It is these white dwarf fluxes which we list in Section 
3.4.21 We then used the white dwarf model atmospheres of 



Gansicke et al.l (| 19951 ) to calculate the white dwarf temper- 



ature T e ff, fixing loggr to the value implied by the zero- 
temperature mass and radius. We then compared these val- 
ues of T e ff an d Mi to the white dwarf cooling models of 
iBergeron etafl (|l995T l in order to find the white dwarf ra- 
dius which is consistent with these values. By comparing this 
to the Eggleton zero-temperature radius, we determined an 
'oversize factor' for the white dwarf, which we found to be 
1.03. Finally, we determine a final white dwarf mass and 
radius using a new mass/radius relation, which is the zero- 
temperature relation scaled by this oversize factor. In theory, 
at this point we should then re-iterate this process and use 
the new mass and radius to refine the temperature determi- 
nation. However, since the first oversize factor is very close to 
1 the effect of further iterations is negligible: the white dwarf 



temperature fit is not affected by the very small change in 
log <j which results from a 3% increase in the white dwarf 
radius. Note also that since the oversize factor is so close 
to 1, while the uncertainty in our temperature measurement 
may be quite large (since it is based on a small number of 
flux measurements) it will not affect our parameter deter- 
minations to a significant degree. 

In the second half of Table [6l we list the binary separa- 
tion a, the masses of the two components, the radial veloc- 
ity semi-amplitudes of the two components and the radius of 
the secondary. This radius is the volume radius of the Roche 
lobe fi lling donor, whi ch we calculate using the approxima- 
tion of lEggletonl (l983h . In our o'-band fits, we find the mass 
of the primary white dwarf to be 0.85 ± O.O4M0. The mass 
of the donor is 0.03 5 ± 0.003Ma and t he donor radius is 
0.047 ± 0.001i? Q . In iMarsh etaU (|2007T I, we reported the 
primary mass to be 0.84 ± O.O5M0 and the donor mass to 
be 0.029 ± O.OO2M0. These values are close to our updated 
findings - the differences are primarily due to the reasons 
discussed in Section 15.11 as well as our accounting for the 



finite temperature of the primary white dwarf in this work. 
Finally, by using our measurements of the white dwarf 



© 0000 RAS, MNRAS 000, 000-000 



16 CM. Copper-wheat et al. 



Table 5. Results from our MCMC fits to the phase-folded light 
curves, as detailed in Section 14,21 q: the mass ratio, i: the bi- 
nary inclination. Hi: the white dwarf radius. Rdisc- the accretion 
disc radius. /3: the power-law exponent for the bright spot. <5: 
the exponent of surface brightness over the accretion disc. I: the 
bright spot scale-length. / c : the fraction of the bright spot taken 
to be equally visible at all phases. <f>: the angle made by the line 
of elements of the bright spot, measured in the direction of bi- 
nary motion from the line tangental to the accretion disc, tfi: the 
angle away from the perpendicular at which the light from the 
bright spot is beamed, measured in the same way as </>. A com- 
pl ete description of all these parameters is given in the appendix 
of ICopperwheat et aljfeoio] 
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contributi on, along with the theo retical absolute magnitudes 
from the iBergeron et al.l (|l995l ) and iHolberg fc Bergeron! 
(2006) cooling models, we can estimate the distance mod- 
ulus for this system. Using the absolute magnitudes for a 
white dwarf temperature of 17, 000K and mass of 0.8M Q , 
our flux measurements imply a distance of 460 - 470pc. 

5.3 AM CVn formation scenarios 

A donor mass of 0.035 ± O.OO3A/0 implies the donor is 
only partially degenerate: a fully degenerate donor in a 
system with this pe riod would have a mass of ~0.020Mq. 
iRoelofs et ail (|2007h measured the masses of five systems us- 
ing parallax measurements obtained with HST, and found 
four of the five to be consistent with a partially degenerate 
donor. The degree of degeneracy varied, but the donor was 
typically found to have a mass of between 2 and 4 times the 
mass of a fully degenerate donor. The lowest estimate was 
for HP Lib, which was found to have a donor mass of 1.6 - 
2.9 times the fully degenerate mass. At 1.75 times the fully 
degenerate mass, SDSS 0926 is by comparison at the lower 
end of these estimates. 

We now examine the finding of a donor mass of 0.035 ± 
0.003Afo in the context of the evolutionary history of this 
system. There are three proposed formation paths for AM 
CVn binaries and all three are consistent with a donor 
that is partially degenerate to some degree. The 'white 
dwarf channel' (|Nelemans et al.l [2001) suggests detached 
close double white dwarfs which are brought into contact 
as a result of angular momentum loss due to gravitational 
wave radiation ( GWR) . INelemans et ail l|200ll ) used a zero- 
temperature donor in their formulation, but iDelove et al.l 
(2005) argued that the donors could be semi-degenerate 
to some degree, depending on the contact time of the bi- 
nary. ThB_se£on^_lOTmation path is the 'helium star chan- 
nel' (llben fc Tutukovlll99ll ). In this scenario the donor is 
a low-mass, non-degenerate helium burning star. Follow- 
ing contact, material is accreted from the helium star onto 
the primary white dwarf until a donor mass of ~ 0.2M© is 
reached, at which point core helium burning ceases and the 
star becomes semi-degenerate. Further mass transfer driven 
by GWR sees the donor mass decrease and the orbital pe- 
riod increase to values consistent with the observed AM CVn 
population. The third scen ario is the 'evolved CV channel' 
(|Podsiadlowski et al.|[2003h which suggests the progenitors 
of AM CVns are CVs with evolved secondaries. The donor 
in this channel is initially non-degenerate and hydrogen-rich, 
but becomes degenerate and helium-rich (but still with a few 
per cent hydrogen) during its evolution before Roche lobe 
over flow. 

INelemans et al.l (|200ll ) approximated the mass-radius 
relationship for a helium star donor and thus modelled the 
evolution of AM CVns formed via the helium star channel. 
Based on the examples provided, an AM CVn with a pe- 
riod equal to that of SDSS 0926 should have a donor mass 
of ^0.05.Mp) if it were for med by this channel. Table 1 of 
Podsiadlo wski et al.l l|2003h lists model parameter determi- 
nations for six AM CVns assuming they were formed via the 
evolved CV channel. For the two systems in this table clos- 
est in period to SDSS 0926 (V803 Cen and CP Eri), these 
results predict the donor mass to be ~O.O4M0 at the point 
at which the orbital period of the system begins to increase. 
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Table 6. Binary parameters for SDSS 0926. We list the mass ratio q, inclination i, and primary radius Ri as determined from the model 
fits. We derive the remaining binary parameters (the binary separation a, the component masses Mi and M2 and the radial velocity 
semi-amplitudes K\ and K2) from these values a nd our cal c ulated mass/radius relation for the white dwarf (Section I5,2[ l, The donor 
radius R2 is calculated using the approximation of lEggletonl jl983T> . 
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Figure 15. Variation in the accretion disc radius over the course 

of our 2006 observations. We divide the three nights of data into 

groups of five or six orbital cycles, which we phase-fold and fit. 

We plot here the accretion disc radius Rdisc scaled by the binary 

separation a, against the MJD of the mid-point of each block of 

data. The three panels of the plot show the results for the first, 

second and third night. 

Since the donor mass in iMarsh et all (|2007f ) was found 
to be 0.029 ± O.OO2M0, lower than the theoretical values 
proposed for the h elium star and evolved CV channels, 
iDelove et all (|2007t) argued that this was evidence for for- 
mation via the white dwarf channel. The updated donor 
mass value we report here of 0.035 ± 0.003Mq is consistent 
with formation via the evolved CV c hannel, and is within 
a few a of the Nel emans et alj (|200ll ) value for the helium 
star channel. We conclude therefore that our current find- 
ings do not strongly preclude any of the three formation 
channels. More precise determinations are likely to be pos- 
sible if we are able to make more observations of SDSS 0926 
in its non-superhumping state, since even after subtracting 
the superhump from our data there is likely to still be some 
residual systematic effect. The discovery of new eclipsing 
AM CVn systems is also key, particularly at the short end 
of the period distribution where there is the biggest discrep- 
ancy between the various donor mass predictions. 

5.4 The disc radius 

In Table[5]we list the apparent accretion disc radius for each 
night as determined from our MCMC fits. This is determined 
from the position of the bright spot, since we assume the 
bright spot to be at the outer edge of the accretion disc. We 



see in this table that the disc radius changes by between 10 
and 20% from night to night. As we discussed in Section|4]2] 
this is also apparent in the lightcurves (Figure [13} with the 
phase of the bright spot eclipse changing between nights. 

The changes we observe here are not due to radial vari- 
ations in a circular disc, rather the disc is non-circular and 
the measured positions of the bright spot sample the possible 
range in disc radii. The superhumps we observe are gener- 
ally taken to imply an elliptical and precessing disc, but the 
disc shape may be more complex than this, with detailed 
numerical simulations suggesting an irregular disc shape i n 
superhumping AM CVn systems (|Simpson fc Woodlll998h . 
The changing bright spot location we find in this system 
is in strong contras t with the study of AM CVn itself by 
iRoelofs et all (2006), who found very little variation in the 
bright spot position. This difference may be related to the 
very different Ms which would be expected in these two sys- 
tems. However, AM CVn was also found to be inconsistent 
with the Patterson e — q relation (Section I3.2.2|> . which is 
presumed to be independent of M. 

We investigated our bright spot findings further by di- 
viding our data into sections, in order to see if the disc radius 
variations can be observed over the course of a single night. 
We split each night of the 2006 data into sections that are 
5/6 orbital cycles in length, and phase- folded and fitted each 
section individually (Figure [15|) . We find there is a notice- 
able increase in radius when we compare the two halves of 
the first night, with R\/a increasing from 0.318 ± 0.007 to 
0.356 ± 0.006. During the second night the disc radius ap- 
pears to be approximately constant over the course of our 
8h observation, and over the course of the third night a 
slight decrease in radius is observed when the two halves of 
the night are compared. These variations appear consistent 
with the precession period of the superhump. 



6 CONCLUSIONS 

In this paper we have presented high time resolution ob- 
servations of the eclipsing P = 28 min AM CVn binary 
SDSS 0926 obtained with the fast CCD camera ULTRA- 
CAM mounted on the William Herschel Telescope. The 
primary aim of these observations was precise parameter 
determinations for the two binary components, using the 
photometric method. We determine the mass ratio to be 
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q = M a /Mi = 0.041 ± 0.002, and the inclination to be 
82.6 ± 0.3 deg. We calculate the mass of the primary white 
dwarf to be 0.85±0.04Mq , and find the donor to be partially 
degenerate with a mass of 0.035 ± O.OO3M0. We also mea- 
sure the eclipse timings with precision, and should be able 
to detect the period change due to gravitational wave losses 
with a third epoch of high time resolution observations. 

We observed in our 2006 WHT/ULTRACAM data the 
superhump that has previously been reported in this source. 
We determine the period of this variation and find it to be in 
agreement with the period excess / mass ratio relationship 
proposed bv lPatterson et al] {2005). This phenomenon is not 
present in the 2009 WHT/ULTRACAM data, but we do see 
evidence for it in some LT/RISE data collected over the 
first half of 2009. In addition to the superhump we observe 
a quasi-periodic oscillation in the 2006 data, with a period 
of around 50s. Another interesting feature of our data is that 
we observe the accretion disc radius to be highly variable. 
We ascribe this to the tidal instability of the outer disc and 
observe changes of up to 10% in the radius over successive 
nights. 

We obtained X-ray and ultraviolet observations with 
XMM-Newton. We found no clear evidence for an eclipse 
in the X-ray light curve. The eclipse is detected in the ul- 
traviolet with the XMM-Newton Optical Monitor, and by 
fitting our model to these data we inferred the ultraviolet 
fluxes of the primary white dwarf. We used these along with 
our ground-based determinations of the white dwarf optical 
colours to determine the temperature of the white dwarf, 
which we found to be 17, 000K. 

Using data collected with the Catalina Real-Time Tran- 
sient Survey, we examined the outbursting behaviour of 
SDSS 0926 over a four and a half year period. We observed 
six outbursts over this period, in which the source flux in- 
creases by ~2 - 2.5 magnitudes. The average time between 
outbursts is ~100 - 200 days. We observed in detail one addi- 
tional outburst with LT/RISE. These data show an increase 
in source flux of greater than a factor of 3, and a rapid decay. 
A lightcurve obtained one day after the detection of the out- 
burst shows a complex structure, with the primary eclipse 
suggesting an asymmetrical distribution of disc flux or possi- 
bly an enhanced bright spot emission. There is additionally 
a dimming prior to the main eclipse which we suggest may 
be due to a warp in the disc itself. 
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